Quasiparticle density of states of rf-wave superconductors in a disordered vortex lattice 
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We calculate the density of states of a disordered inhomogeneous d-wave superconductor in a 
magnetic field. The field-induced vortices are assumed to be pinned at random positions and the 
effects of the scattering of the quasi-particles off the vortices are taken into account using the singular 
gauge transformation of Franz and Tesanovic. We find two regimes for the density of states: at very 
low energies the density of states follows a law p(e) ~ po + P\e\ a where the exponent is close to f . 
A good fit of the density of states is obtained at higher energies, excluding a narrow region around 
the origin, with a similar power law energy dependence but with a close to 2. Both at low and at 
higher energies po scales with the inverse of the magnetic length ~ \fB)- 

PACS numbers: 74.25.Qt, 74.72-h 



The effect of disorder on the low-energy density of 
states of d~wave superconductors has been a subject of 
considerable recent interest both for practical and for the- 
oretical reasons Q . From the practical point of view the 
presence of disorder pinning mechanisms is important to 
prevent energy dissipation due to the motion of the vor- 
tices in an external field. From the theoretical point of 
view several conflicting predictions have appeared in the 
literature. Some progress toward understanding the dis- 
parity of theoretical results has been achieved realising 
that the details of the type of disorder affect significantly 
the density of states [lj. In contrast to conventional 
gapped s-wave superconductors, the presence of gapless 
nodes in d-wave superconductors leads to power law be- 
havior for the low-T thermodynamic properties and is 
expected to affect the transport properties. Linearizing 
the spectrum around the four Dirac-like nodes it has been 
suggested that the system is critical. It was obtained that 
the density of states is of the type p(e) ~ |e| a , where a is 
a non-universal exponent dependent on the disorder, and 
that the low energy modes are extended states (critical 
metal) [j| . Taking into account the effects of inter-nodal 
scattering (hard-scattering) it has been shown that an 
insulating state is obtained instead, where the density of 
states still vanishes at low energy but with an exponent 
a = I independent of disorder |3| . The addition of time- 
reversal breaking creates two new classes designated spin 
quantum Hall effect I and II, due to their similarities 
to the usual quantum Hall effect, corresponding to the 
hard and soft scattering cases, respectively The pro- 
posed formation of a pairing with a symmetry of the type 
d + id breaks time-inversion symmetry Q but up to now 
remains a theoretical possibility. On the other hand ap- 
plying an external magnetic field naturally breaks time- 
reversal invariance and therefore it is important to study 
the density of states in this case. 

Moreover, the interaction between the superconductor 
quasiparticles and the vortices induced by the external 
magnetic field has also been a subject of considerable 



debate 0, IS • I n the presence of the vortices the quasi- 
particles feel the combined effect of the external mag- 
netic field and of the spatially varying field of the chiral 
supercurrents. Performing a gauge transformation to ef- 
fectively reduce the system to one in a zero average mag- 
netic field it was shown Q that the natural low energy 
modes are Bloch waves rather than the Dirac Landau lev- 
els proposed in U @ . These results also showed that the 
quasiparticles besides feeling a Doppler shift caused by 
the moving supercurrents ||| (scalar potential) also feel 
a quantum "Berry" like term due to an Aharonov-Bohm 
scattering of the quasiparticles by the vortices (vector 
potential). 

In this Letter we will be concerned with the effect of 
positional vortex disorder on the density of states of a su- 
perconductor in an external magnetic field. The density 
of states of a dirty but homogeneous s-wave superconduc- 
tor in a high magnetic field where the scattering of the 
quasiparticles off scalar impurities was considered using 
a Landau level basis jfj. For small amounts of disor- 
der it was found that p{e) ~ e 2 but when the disorder is 
higher than some critical value a finite density of states is 
created at the Fermi surface. In the same regime of high 
magnetic fields but with randomly pinned vortices and no 
impurities the density of states at low energies increases 
significantly with respect to the lattice case suggesting a 
finite value at zero energy [ljj. Refs. considered 
the effects of random and statistically independent scalar 
and vector potentials on d-wave quasiparticles and it was 
predicted that at low energies p(e) ~ po + ae 2 , where 
po ~ B 1 / 2 . To our knowledge the effect of randomly 
pinned discrete vortices on the spectrum of a d-wave su- 
perconductor has not been addressed previously. 

We will consider the lattice formulation of a disor- 
dered d-wave superconductor in a magnetic field. We 
start from the Bogoliubov-de Gennes (BdG) equations 
Hip — ei/j where = (u*(r),v*(r)) and where the 
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matrix Hamiltonian is given by 

«-(it.i) 

with j3, E3 

A = Ao^e^M^e^W. (1) 

The sums are over nearest neighbors (8 = ±x, ±y on the 
square lattice); A(r) is the vector potential, $su(r) = 
u(r + S), fjs = 1/4 f° r s-wave pairing and i)$ = (—l) Sv §s 
for d-wave pairing. In eq. (1) we have factorized the 
phase of the order parameter and have taken the Lon- 
don limit assuming that the amplitude Ao is constant 
everywhere in space, which is valid in the regime of low 
fields where the size of the vortex cores is negligible. It is 



convenient to perform a singular gauge transformation to 
eliminate the phase of the off-diagonal term in the matrix 
Hamiltonian. In such a way the magnetic field is com- 
pensated by an array of opposing half-fluxes. We carry 
out the unitary transformation |jj H. — > U~ l 7iU , where 

/ e i<Mr) o \ 
U ~\ e-***M J 

with </>a(i") + 0e(r) = <X r )- The vortices are separated 
into two groups A and B positioned at {Yf}i=\.N A and 
{rf }i=i.N B such that the vortices A are only visible to 
the particles and the vortices B are only visible to the 
holes. The phase fields 4>^=a,b are defined through V x 
V<^ p (r) = 2nzJ2id( r — r f)- The resulting system is 
effectively in a magnetic field = B— cf) z ^ ■ 5(r— rf) 
that vanishes on average if Na — Nb ■ Under the unitary 
transformation the BdG equations convert in the d-wave 
case to [l3( 



-<^e lV ^ (r) w(r + <5) -e F u(r) + A ^ e l ^ (r) (-l)^(r + (5) = eu(r) 

8 8 

A ^ e ,ii(r) (-l) J »ii(r + ^) + ^ e- !V ' (r) i)(r + + £ F i)(r) = ev(r) (2) 

I 



where the phase factors are given by [ijjf V^(r) = 
J r r+6 ■ dl and As(r) = \ J^ +S (kf - kf ) • dl. Here 
hk^ = ?7iv^ = h~V(j)^ — |A, where k^ is the superfluid 
wave vector for the ^t-supercurrent. This quantity can be 
calculated for an arbitrary configuration of vortices 0] 
like 

i— 1 

Here A is the magnetic penetration length and the sum 
extends over all vortex positions. The BdG equations are 
then solved taking an arbitrary configuration of vortices 
and the density of states is calculated in the standard 
way. The average over disorder is then carried out calcu- 
lating the density of states for each vortex configuration 
and then performing the average on the density of states 
over 100 different configurations. 

The situation where the vortices are regularly dis- 
tributed in a lattice was treated before 0. Since the 
average magnetic field vanishes it is possible to solve the 
BdG equations using a standard Bloch basis since the 
supercurrent velocities are periodic in space and there is 
no need to consider the magnetic Brillouin zone. Taking 
the continuum limit and linearizing the spectrum around 
each node effectively decouples the nodes. It was shown 



that the low-energy quasiparticlcs arc then naturally de- 
scribed as Bloch waves and not Dirac-Landau levels as 
previously proposed 0, 0]- However, it was found that 
in the linearized problem different assignments of the A 
and B vortices lead to somewhat different spectra, which 
was unexpected |13(. It was found that taking the the- 
ory on the lattice regularized this problem and indeed 
the system has a manifest internal gauge symmetry such 
that the spectrum is independent of the A-B vortex as- 
signments, as it should be. Moreover, the lattice formu- 
lation explicitly involves internodal contributions which, 
as discussed above, are important for the properties of 
the density of states in the disordered case. In the vortex 
lattice case, however, it was found that only in special 
commensurate cases (for the square lattice) the inclusion 
of the internodal contributions is relevant since only in 
such cases a gap develops due to the interference terms 
between the various nodes, estimated to be of the order of 
VB. In a general incommensurate case the interference is 
not relevant leading to qualitatively similar spectra. The 
spectrum is gapless with a linear density of states at low 
energy |l3| . One would therefore expect that in a general 
disordered vortex case internodal scattering might not 
be relevant (particularly for high Dirac cone anisotropy 
t/A ). 
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FIG. 1: Quasiparticle density of states for different magnetic 
fields B = 1/2 (dashed line) and B = 1/200 (solid line) in 
units of hc/(2e8 2 ). The inset shows the Landau level quanti- 
zation at high-energy for B — 1/50. 
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FIG. 2: Quasiparticle density of states for different magnetic 
fields B = 1/200 (A), B = 3/200 (A), B = 5/200 (■), B = 
7/200 (□), B = 9/200 (•), B = 11/200 (o), B = 20/200 (♦) 
and B = 25/200 (0) in units of hc/{2eS 2 ). The linear system 
size is L — 20 and the parameters are /i = — 2.2i, A<> = t. 
The inset shows the fits of the density of states (solid lines) 
inside the presented energy interval. For clarity not all the 
fits are shown. 



We consider a square lattice with lattice constant 5 
(taken as 1), where the electron hopping is described by 
Hamiltonian (1). Application of the external magnetic 
field generates vortices which are placed at the center of 
a plaquette (unit cell). The number of vortices in the 
two-dimensional system of size L x L is proportional to 
the quantized magnetic flux piercing through the system. 
We parametrize the intensity of the magnetic field by 
the ratio of the number of vortices, N^, (divided equally 
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FIG. 3: Density of states po near the Fermi energy and expo- 
nent a (Inset) extracted from fits for different linear system 
sizes L = 10 (A), L = 12 (o), L = 16 (□), and L — 20 (o). 
The fit of po gives a square root dependence in B shown by 
the solid line. 



in two groups A and B) by the number of lattice sites 
B = L ^ L . The vortices are distributed randomly over 
the L x L plaquettes. We consider then periodic boundary 
conditions and solve the BdG equations numerically. We 
consider A — > oo. In this limit the repulsive interaction 
between vortices is not screened and therefore the vortex 
distribution is not strictly arbitrary. We assume that the 
pinning centers are strong enough to overcome the vortex 
repulsion. We have checked that the case with a random 
distribution of vortices is qualitatively the same as for the 
case where the vortex positions are allowed to vary with 
a radius of a few unit cells around a lattice position. We 
find that, as the disorder is introduced, the low energy 
density of states generally increases. 

In Figure 1 we show the density of states for the cases 
of a weak (B = 1/200) and of a strong magnetic field 
(B = 1/2; for this high value the density of vortices is 
high and strictly we are in a regime where the size of the 
vortex cores can not be neglected) . In the case of a strong 
magnetic field the Landau level structure at high energies 
is clear and it extends to low energies superimposed by 
the effects of disorder and the effect of the low-energy 
modes close to the <i-wave nodes. At weak fields the 
density of states is small at low energies having a dip 
close to zero energy. We have checked for finite size effects 
on the spectrum. For system sizes larger than 16 x 16 
the density of states at not very low energies converges 
and the finite size dependence is negligible. In Figure 
2 we show the density of states for a system with size 
20 x 20 and for various magnetic fields. The density of 
states at small energies is finite up to quite small energies 
where there is a dip to a value that decreases as the 
magnetic field decreases. Only for quite small magnetic 
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FIG. 4: Low-energy quasiparticle density of states p(e) for 
different magnetic field B = 1/200 (A), B = 3/200 (A), B = 
5/200 (■), B = 7/200 (□), B = 9/200 (•), B = 11/200 (o), 
B = 20/200 (♦) and B = 25/200 (0) in unit of hc/(2e6 2 ). 
For each data sets the linear system size is L — 20, the same 
as in Fig. [5] In the left panel the solid lines are linear fits of 
the dip region below e = 0.02t of the type p(e) = podip +/3|e|. 
In the right panel we present the near scaling at low energies 
(see text). 

fields the density of states approaches zero at the origin. 
Neglecting the narrow region close to the origin we have 
fitted the density of states using the power law 

p(e)=po + P\e\ a . (4) 

In the inset of Figure 2 we show the fits for the various 
values of the magnetic field. Reasonable fits are obtained 
taking a ~ 2 and we obtain that po ~ B 1 / 2 . In Figure 3 
we show the magnetic field dependence of the parameters 
of the fit for different system sizes L = 10, 12, 16, 20. The 
various system sizes fit in the same universal curve indi- 
cating that the finite size effects are negligible. Note that 
in the lattice case the density of states at low energies is 
linear 0, El E-l E3| (this result differs from the behavior 
obtained by others for a d-wave superconductor with no 
disorder H,[lg, where p(e ~ 0) ~ B 1 / 2 ). The finite den- 
sity of states at zero energy is therefore a consequence of 
finite disorder. 

In Figure 4 we focus on the narrow region close to 
e = for the same set of parameters considered in Fig- 
ure 2. Except for the lowest field case B = 1/200 the 
density of states seems to be finite at zero energy. The 
particular field density of B = 1/200 corresponds to only 
two vortices (since the size of the system is 20 x 20). 
As shown in [l4| in this case the spectrum is usually 
gapped and therefore the density of states vanishes at 
zero energy. Performing a fit like in cq. (4) we obtain 
an exponent which is now close to 1. In this regime 
Podip also scales with \J~B and the slope scales linearly 



FIG. 5: Slope j3 and density of states podip (Inset) close to 
e ~ extracted in the dip region from data with different 
linear lattice sizes L = 10 (A), L = 12 (♦), L = 16 (■) and 
L — 20 (•). Dashed lines are linear fits of /3. In the inset podip 
is shown for L — 20 (o). The fit of podip gives a square root 
dependence in B shown by the solid line. 



with B. In this low energy regime the finite size ef- 
fects are still noticeable but the dependence on the mag- 
netic field is common to the various system sizes. At 
these low energies the density of states for the vari- 
ous system sizes appears to be of the following approxi- 
mate form p{e) ~ {1/uj h ){1/1 2 )T ((e/Lu H )(5 2 /I 2 )) , where 
u>h ~ \/AqB and J 7 is a universal function. In the left 
panel of Figure 4 we show p(e) for various fields while 
in the right panel we illustrate the near scaling at low 
energies consistent with T{x) ~ c\ + cix at small x. 
In Figure 5 we show the field dependence of the slope 
of p(e) and of the zero energy density of states. Note 
that f3(B — > 0) (Fig. 5) appears to be small but fi- 
nite, consistent with a crossover to a "Dirac node" scaling 
(p(e) ~ (l/w H )(l/P)J r (e/w_ ff )) at very low fields. 

In summary, we have calculated the density of states 
of a disordered d-wave superconductor in a pinned vortex 
array. Both the disorder and the magnetic field fill the 
density of states at low energies. In general we find a fi- 
nite density of states at zero energy except in the limit of 
very small magnetic fields. The density of states deviates 
from the zero energy value by a power law. Excluding a 
narrow region close to zero energy a good fit is obtained 
with an exponent close to 2. Performing the same type 
of fit at very low energies a good fit is obtained with an 
exponent close to 1. In the very low magnetic field limit 
the dip of the density of states is more pronounced. The 
zero energy density of states scales with the inverse of the 
magnetic length (V~B). Except for the zero energy finite 
value the energy dependence of the density of states in 
the case with disorder is similar to the lattice case. This 
suggests that the vortex disorder does not dramatically 
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affect the density of states at low energies. A prelimi- 
nary analysis of the inverse participation ratio indicates 
that the states are still extended, as in the lattice case. 
Our results for the lowest and zero energies are clearly 
most susceptible to finite size effects- we will report on 
a detailed study of this region separately. In the gapped 
s-wave case however the disorder introduces states in the 
gap thereby changing qualitatively the low energy den- 
sity of states, as in the high field limit [inj . These results 
will be presented elsewhere. 
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